NSWCCD-70-TR-2000/038  Structural  Fuzzies  At  Sea 


Carderock  Division 

Naval  Surface  Warfare  Center 

9500  MacArthurBlvd. 

West  Bethesda,  MD  20817-5700 


NSWCCD-70-TR-2000/038  February  2000 

Signatures  Directorate 

Research  and  Development  Report 


Structural  Fuzzies  At  Sea 

by 

G.  Maidanik 


I 


Approved  for  public  release;  Distribution  is  unlimited 


*«•«»,.  20000509  198 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 


REPORT  DOCUMENTATION  PAGE 


la.  REPORT  SECURITY  CLASSIFICATION 
Unclassified 


2a.  SECURITY  CLASSIFICATION  AUTHORITY 


2b.  DECLASSIFICATION/DOWNGRADING  SCHEDULE 


Form  Approved 
OMB  No.  0704-0188 


1b.  RESTRICTIVE  MARKINGS 
None 


3.  DISTRIBUTION /AVAILABILITY  OF  REPORT 

Approved  for  public  release. 

Distribution  is  unlimited. 


4.  PERFORMING  ORGANIZATION  REPORT  NUMBER(S) 
N  S  WCC  D-70  -TR — 200  0/0 38 


5.  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 


6a.  NAME  OF  PERFORMING  ORGANIZATION  6b.  OFFICE  SYMBOL  7a.  NAME  OF  MONITORING  ORGANIZATION 

(if  applicable) 

Carderock  Division,  NSWC  Code  7030 


6c.  ADDRESS  (City,  State,  and  ZIP  Code) 

9500  MacArthur  Blvd. 

West  Beth esda,  MD  20817-5700 


8a.  NAME  OF  FUNDING/SPONSORING 
ORGANIZATION 


8b.  OFFICE  SYMBOL 
(if  applicable) 


7b.  ADDRESS  (City,  State,  and  ZIP  Code) 


9.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 


8c.  ADDRESS  (City,  State,  and  ZIP  Code) 


11 .  TITLE  (Include  Security  Classification) 
Structural  Fuzzies  At  Sea 


12.  PERSONAL  AUTHOR(S) 
G.  Maidanik 


13a.  TYPE  OF  REPORT 
Research  &  Development 


13b.  TIME  COVERED 

FROM  010300  TO  030100 


16.  SUPPLEMENTARY  NOTATION 


10.  SOURCE  OF  FUNDING  NUMBERS 


PROGRAM 
ELEMENT  NO. 


PROJECT 

TASK 

NO. 

NO. 

WORK  UNIT 
ACCESSION  NO. 


14.  DATE  OF  REPORT  (Year,  Month,  Day)  15.  PAGE  COUNT 
2000  March  1  54 


COSATI  CODES 


GROUP 


SUB-GROUP 


18.  SUBJECT  TERMS  (Continue  on  reverse  if  necessary  and  identify  by  block  number) 

Statistical  Energy  Analyses  (SEA) 

Structural  Fuzzy 
Modal  Density 
Structural  Damping 


18.  ABSTRACT  (Continue  on  reverse  if  necessary  and  identify  by  block  number) 

A  rudimentary  statistical  energy  analysis  (SEA)  can  be  employed  to  decipher  the  functioning  of  structural  fuzzies.  A  structural  fuzzy  is  a  device  that,  when 
coupled  strongly  to  a  structure,  increases  the  damping  of  that  structure.  A  structural  fuzzy  must  be  endowed  with  a  modal  density  and  structural  damping 
that  by  far  exceeds  those  of  the  structure  to  which  it  is  strongly  coupled. 


20.  DISTRIBUTION/ AVAILABILITY  OF  ABSTRACT  21.  ABSTRACT  SECURITY  CLASSIFICATION 

□  unclassified/unlimited  □  sameasrpt.  □  OTIC  USERS  Unclassified 


[  22a.  NAME  OF  RESPONSIBLE  INDIVIDUAL 
Maidanik,  G. 


JD  Form  1473,  JUN  86 


22b.  TELEPHONE  (Including  Area  Code) 
(301)  227-1292 


Previous  editions  are  obsolete 


22c.  OFFICE  SYMBOL 
NSWCCD  7030 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 


ABSTRACT 


A  rudimentary  statistical  energy  analysis  (SEA)  can  be  employed 
to  decipher  the  functioning  of  structural  fuzzies.  A  structural  fuzzy  is  a 
device  that,  when  coupled  strongly  to  a  structure,  increases  the  damping  of 
that  structure.  A  structural  fuzzy  must  be  endowed  with  a  modal  density 
and  structural  damping  that  by  far  exceeds  those  of  the  structure  to  which 
it  is  strongly  coupled. 

Friday,  October  15,  1999  1 1:00 
Room  210 
Photonic  Center 
8  St.  Mary’s  Street 
Boston,  MA  02215 

This  report  is  based  on  the  seminar  given  at  Boston  University  as 
stated  above.  Viewgraphs  are  presented  on  the  right  and  comments 
relevant  to  that  viewgraph  are  placed  on  the  left.  A  few  of  the  remarks 
made  by  the  audience  are  included  without  reference.  Nonetheless,  the 
author  is  thankful  for  these  participations  by  the  attendees. 


l 


A  Note: 


Recently  a  full  fledged  effort  has  been  launched  to  develop  analyses  which  yield 
information  transcending  that  which  is  issued  by  SEA.  In  particular,  there  are  claims  that 
the  information  issued  by  SEA  in  the  low  and  mid-frequency  ranges  is  deemed 
inadequate  for  certain  response  descriptions  and  noise  control  purposes.  These  issues 
were  deliberately  forfeited  by  assumptions  imposed  during  the  development  of  SEA  and, 
therefore,  they  lie,  by  definition,  beyond  the  scope  of  SEA.  The  assumptions,  so 
conceived,  render  SEA  a  more  tractable  analysis.  It  is,  therefore,  impertinent  to  blame 
SEA  for  any  inadequacies  resulting  from  these  forfeits.  Moreover,  there  are  those  who 
interpret  the  results  in  the  low  and  mid-frequency  ranges  by  the  failure  of  SEA  to  yield 
information  that  is  suppressed  apriori  in  the  development  of  SEA.  Yet,  by  definition, 
SEA  can  and  will  yield  the  information  that  it  is  designed  to  yield,  independent  of  the 
frequency  range.  Notwithstanding  that  in  a  given  frequency  region,  the  deviation  from 
the  mean-value  may  fluctuate  more  and  more  as  the  violation  of  the  condition  of  the 
modal  overlap  becomes  more  severe  [1,2].  In  SEA  the  mean-value  is  the  fundamental 
description.  Indeed,  the  very  strength  of  SEA  is  that  it  chooses  the  right  thing  (details)  to 
be  ignorant  about  as  a  way  to  learn  something  (mean-values).  It  takes  chutzpah  to 
ascertain  the  in  the  low  and  mid-frequency  ranges  by  the  apparent  failure  of  SEA  to 
account  for  phenomena  that  lie  outside  its  scope.  It  transpires  that  the  low  and  mid¬ 
frequency  ranges  are  often  determined  by  the  failure  of  those  who  profess  to  understand 
what  SEA  says  is,  is ! 
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Revisiting  Elements  in  Statistical 
Energy  Analysis  (SEA) 

G.  Maidanik,  NSWCCD  (DTMB) 

AGENDA 

•  A  brief  review  of  SEA  is  presented. 

•  The  modifications  in  SEA  due  to  the  influence  of  strong 
couplings  among  the  substructures  are  discussed. 

•  Some  remedial  modifications  to  SEA  are  heuristically 
proposed. 

•  Two  noise  control  problems  are  described  and  interpreted 
in  terms  of  SEA: 

1.  Control  of  loss  factors  (dissipations)  in  pipe-like 
structures  containing  beads.  The  beads  are  modeled 
to  be  a  fluid  of  low  sound  speed  and  of  high  damping. 

2.  Control  of  vibrations  by  increasing  damping  and 
isolation.  In  particular,  a  note  is  taken  of  the 
supplementarity  of  damping  and  isolation  in  vibration 
control. 
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Mathematical  Symbols 
Scalar 

Vector  of  Rank  N 


v 


Transposed  v+  =  {vl5  v2, _ v-. . .  .v^} 

Complex  conjugate  v*  =  {v*} 

v+  v*  =  a  Scalar  ;  ^  v  =  a  Matrix 


Matrix  of 
Rank  V 


(*11  *12 
*21  *22 


[*N1 


e.g., 

Impedance  Matrix  Eq.  Admittance  Matrix  Eq. 
Z  v  —  F  v  =  y  F 

where  symbolically  ^  =  (Z)_1 


V.2 
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Viewgraph  3 


The  division  of  a  complex  structure  into  constituent  substructures  is  crucial  in 
SEA.  Yet,  rarely  is  this  division  discussed,  and  specifications  for  it  are  rarely  laid  down. 
In  part,  the  reluctance  stems  from  the  belief  that  “An  arbitrariness  (non-uniqueness)  in  the 
service  of  a  good  approximation  is  not  a  vice,  it  is  a  useful  device”  and,  in  part,  apriori 
rules  are  difficult  to  state.  In  addition,  SEA  is  a  function  of  the  frequency  (steady  state) 
variable  (CO)  or  the  temporal  (transient  state  variable  ( t ) .  The  division  may,  therefore, 
be  dependent  on  ( co )  or  ( t ) ;  e.g.,  a  division  at  (COl )  may  be  different  from  a  division  at 
( co2  )  where  COl  ^  C02-  That  dependence  of  the  division  is  seldom  discussed  but  its 
significance  is  here  recognized.  In  the  final  analysis,  it  turns  out  that  experience  is  often 
the  best  guide  to  this  division.  Needless  to  state,  the  larger  the  number  of  substructures 
the  higher  is  the  frequency  at  which  the  fundamental  resonance  frequency  lies,  and  the 
higher  is  the  rank  of  the  matrix  equation  of  motion  at  SEA.  Therefore,  the  higher  the 
number  of  substructures  the  higher  is  the  frequency  of  validity,  and  the  more  difficult  is 
the  manipulations  of  the  equation  of  motion.  On  the  other  hand,  the  smaller  the  number 
of  substructures  the  more  details  are  suppressed;  excessive  suppression  of  details  may 
render  the  equation  of  motion  moot  and  without  substance. 

None  of  these  extremes  are  happy  circumstances.  It  is  hoped,  in  this  connection, 
that  some  idea  for  the  appropriate  division  procedure  will  emerge  in  the  reading  of  this 
report. 


Statistical  Energy  Analysis  (SEA) 


(.Zaa(Xa)Say  jdXyZay(Xa  IXy)(l  3(yy  ^ ) 

% 

is  a  vector  that  spans  the  dimensionality  of  the  (a)th 
substructure.  Thus,  if  the  (cc)lh  substructure  is  one¬ 
dimensional,  is  a  scalar.  If  it  is  two-dimensional,  xa  is  a 
two-vector  and  if  it  is  three  dimensional,  xa  is  a  three-vector. 
The  dependence  of  quantities  on  the  time  variable  (t)  or  on 
the  frequency  variable  (i co )  is  implicit  in  this  scalar  or  matrix 
impedance  equation  of  motion. 


V.3 


5 


Viewgraph  4 

A  generic  division  of  a  complex  structure  into  substructures  is  represented  in 
Viewgraph  (3).  The  choice  of  the  division  must,  however,  be  compatible  with  the 
requirement  that  the  self-impedance  operator  of  the  substructures;  e.g.,  Zaa(xa)  of  the 
(a)  th  substructure,  be  eigenoperators  so  that  a  complete  set  of  orthogonal  mode  shape 
functions  may  be  assigned  to  each  substructure;  e.g.,  (j)aj  (Xa )  are  assigned  to  the  a  th 
substructure  as  depicted  in  Viewgraph  (3).  In  this  manner  the  individual  mode  shape 
vectors  and  the  overall  mode  shape  vector  may  be  defined;  (j) aj  (Xa )  =  {fyaj (Xa)}  and 
<P(X)  =  {0a  (*«)},  respectively. 

The  choice  for  the  self-impedance  operators  and,  therefore,  for  the  coupling 
impedance  operators  and  the  mode  shape  functions  is  not  unique;  they  may  be  chosen  in 
a  variety  of  ways,  however,  each  way  may  be  related  to  the  others.  With  experience  the 
better  choices  may  be  identified. 
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A  complex  composed  of  a  number  of  coupled  dynamic 
systems 

Z(x)  V(x)  =  p  (X) 

£  =  k) 

£(£)  =  kk)} 

&(*)  =  \PeJxa)\ 

^(£)  ^Zaa (Xa)^ay  JdXyZay(xa  lx^)(l  ^^7 ) ) 

The  rank  of  the  equation  of  motion  is  equal  to  (N) ,  the 
number  of  dynamic  systems  in  the  complex. 

The  self-impedances;  e.g.,  Zaa(xa),  are  eigenoperators  in  the 
sense  that 

Z>aa  (-X'a  ^^Paj  (-^a  )  ^aj aj  aj  (-^a  )  ’ 

2«k)  =  ;  ?(2)  =  {$?«(*«)} 


Viewgraph  5 


Employing  the  chosen  complete  and  orthogonal  sets  of  mode  shape  functions,  the 
modal  equation  of  motion  may  be  derived  as  stated  in  Viewgraph  (5).  The  rank  of  this 
modal  equation  may  be  quite  large,  especially  at  the  higher  frequency  ranges. 
Nonetheless,  if  this  equation  can  be  solved  for  the  modal  response  ( v )  and  the  spatial 
dependence  of  the  mode  shape  functions  is  retained,  the  response  Va(Xa),  of  typically 
the  oc  th  substructure,  may  be  restored  and  stated  as  prescribed  in  Viewgraph  (6). 
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Modal  decomposition  — ►  removal  of  the  spatial 
dependence  in  the  equation  of  motion  — ►  modal  equation 
of  motion 

~  JZe  9  ~  ~ccy  —^7^—  ^LacS')  ’ 

£=  tU  ;  ;  Pe  =  {ZeJ  ;  £««  =  ttW  ; 

£aa  (Zqjaj  ^ jl )  ’  ~«7  ^ajyq )  » 

Zqjyq  j dXy  j dxa  aj (xa )Zay (,Xa  I X7 )  tyyqi^Xy) 

The  rank  of  the  modal  equation  of  motion  is  equal  to 
X  AL,  where  N„  is  the  number  of  all  relevant  modes  in  the 

a  u  a 

(a)th  dynamic  system. 
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Viewgraph  6 

Again,  the  modal  equations  of  motion  are  stated  in  Viewgraph  (6).  Since  a  set  of 
mode  shapes  is  assumed  to  be  complete  and  orthogonal,  the  modes  in  each  substructure 
are  not  coupled,  couplings  occur  only  between  modes  belonging  to  different  substructures 
(dynamic  systems). 

The  inverse  equation  of  motion  is  stated  in  order  to  indicate  that  if  (Z )  is  known, 
its  inverse  (y )  may  be  used  to  directly  determine  the  response  vector  (v)  once  the 
external  drive  vector  (p  )  is  given. 

The  high  rank  of  the  impedance  matrix  or  the  admittance  matrix  modal  equation 
of  motion  may  render  a  solution  well  nigh  impossible.  Means  to  reducing  the  rank  may 
thus  be  advantageous,  even  at  the  expense  of  forfeiting  details. 
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Recovery  of  the  spatial  dependence 

Va(xa)  ~  Y,a  *  *£a(Xcc)  ’  Pea(X )  —  Pea *  5?aC*a) 

Thus,  if  the  eigenvectors  <Pa(xa)  are  appropriately  stored, 
they  can  be  used  in  this  reconstitution  when  desired. 

Modal  equation  of  motion 

Z  V  —  P e  \  Z  —  (^aa^cr/-  ^or/O-- 

“  (Zajaftjd  kay  = 

No  modal  couplings  in  Modal  couplings  across 

a  dynamic  system  dynamic  systems 

I  =  J  £  ;  v  =  fc.)  ;  2=<z)-1  :  ft  =  (W 


V.6 
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Viewgraph  7 

The  forfeiting  of  details  starts  with  the  abandonment  of  spatial  information  within 
each  substructure;  only  the  response  of  each  mode  is  registered  without  heed  to  its  spatial 
distribution  -  the  mode  is,  then,  reduced  to  be  a  simple  harmonic  oscillator.  Once  the 
spatial  information  is  forfeited,  the  phase  of  each  harmonic  oscillator  loses  its 
significance.  If  forfeiting  the  phases  can  be  advantaged  in  the  effort  to  reduce  the  rank  of 
the  modal  equation  of  motion,  let  it  be  done.  The  phases  are  forfeited  by  quadratizing  the 
modal  equation,  as  stated  in  Viewgraph  (7).  Then,  the  real  part  of  the  resulting  equation 
is  executed,  again,  as  stated  in  Viewgraph  (7).  The  fundamental  statement  of  SEA  is 
presented.  This  statement  is  extrapolated  from  findings  pertaining  to  the  interaction 
between  two  linear  harmonic  oscillators  [2,3], 
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Return  to  the  modal  equation  of  motion 

5.  ^  ~  ^  e  ’  ^ocjaj^aj  2  ^h^ajyq^yq^  &  ay ) —  P  eaj 

.  7  q 

In  this  equation  <P  (x)is  forfeited  by  spatial  averaging. 
Now  the  phases  in  v  are  forfeited  by  quadratizing  the 
modal  equation 

Zqjqj^qj^  ~~  2  X  (To/)  ^ajyq^yq  0-  ~8ay )  “  O^cgf )  P eaj  * 

y  # 


Take  the  real  part  of  and  normalize 

ReE£  (v<9-)*(Za,-w/fi>)(vw)(l-5ay)}=(?re(9-/G)) 

J  J  \J  q  J 

V" 

+  2  lLPajyq(-P'aj~P'yq)  Q~8ay) 

r  q 


This  is  the  fundamental  statement  of  SEA  with 


Small  if  coaj^coyq 
Large  if  (oaj-(Dyq 


V.7 
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Viewgraph  8 

The  definition  of  the  modal  energy  (£„•)  of  the  ( /’)th  mode  (harmonic 
oscillator)  in  the  (a)th  substructure  is  stated  in  terms  of  the  mass  (Ma)  of  that 
substructure.  The  quantity  (£ Kaj)  is  the  kinetic  modal  energy  stored  in  the  referenced 
mode.  The  real  part  of  {  (V™  )  peaj  }  is  equated  to  the  external  power  input  into  the 
(j)  th  mode  in  the  (a)  th  substructure. 

The  span  function  is  defined  in  terms  of  a  frequency  bandwidth  (A&>)  centered 
about  ( CO )  that  includes  all  modes  within  it.  Thus,  (•  •  •  -)^(o  represents  frequency 
averaging  over  the  quantities  that  pertain  to  the  resonance  frequencies  that  lie  within 
(A (0) .  In  transiting  from  Equation  (a)  to  (b)  on  Viewgraph  (8),  equipartion  of  modal 
energy  is  assumed;  namely,  £aj  -  £aj  =  £a  and  £yq  —  £yf  —  £y  of  modes 
that  lie  within  the  modal  span  functions  U aj  and  Uyq,  respectively.  Equation  (c)  on 
Viewgraph  (8),  merely  generalized  Eq.  (b)  and  cast  this  generalization  in  matrix  form. 
Note  that  ffya  =  ffay  but  T]ya  is  not  necessarily  equal  to  T]ay.  Indeed,  the  two 
quantities;  namely,  r\ya  and  rjay,  are  related  by  a  consistency  relationship,  as  stated  in 
Viewgraph  (9). 


14 


Modal  equation  of  SEA 


^ aj  —  I  Vaj  I  —  2£Kaj 

naj  =  Re{(vajr  peaj\ 

^ lajaj  {~q/q/  ^  I 


Using  the  modal  span  functions  Uaj((D,(Daj,riaj) 

(tlajaj^aj  X  X ^7a/yg (^q/'  ~ £yq ^0-~Say)  —  (7Teaj  / CO)) ^  ,  (a) 


^lacc^a  X  Vya^a  £y)0-  5ay)  —  (jCeal (O)  ,  (b') 

7 


(; ne ; 

®  =  (£0<5ar)  ; 


^  “  (^a  &ay  ~  Vya  0-  ~  <5ay  ) )  . 

€  =  {£«}  J  &  =  foal 

^ya  —  ^yVya  >  Hya  —  (.Vytjaj')  Aco 


V.8 
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Viewgraph  9 


Global  density  quantities  are  defined  in  terms  of  the  (typical)  modal  quantities  and 
the  modal  densities.  The  modal  density  na  ((d)  of  the  (a)  th  substructure  counts  the 
modes  per  unit  frequency  that  reside  about  the  frequency  ( CO ).  Thus,  na(co)(Aco)  are 
the  number  of  modes  that  are  spanned  within  the  modal  span  function.  This  number  is 
(Uaj (CO, COaj , riaj)Aco .  [cf.  Viewgraph  (8).]  As  indicated  on  the  figure,  a  substructure 
at  SEA  is  completely  defined  in  terms  of  its  loss  factor;  e.g.,  (T\aa )  for  the  (a)  th 
substructure,  of  its  modal  density;  e.g.,  ( na )  for  the  (a)  th  substructure,  of  its  mass; 
e.g.,  ( Ma )  for  the  (a)th  substructure*  and,  finally,  of  its  coupling  loss  factors;  e.g., 

T]ay  for  the  coupling  between  the  (7)th  and  the  (a)th  substructures.  Consistent 
relationships  exist  among  the  coupling  loss  factors  as  stated  in  Viewgraph  (9). 


*  The  mass  of  a  substructure,  (Ma )  for  the  (a)  th  substructure,  is  seldom  stated 
explicitly,  but  is  always  implied. 
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(Vyal^ay)  =  (nr  / na)  Consistency  Relationship 


I 

I 

I 

I 


V.9 
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Viewgraph  10 


In  Viewgraph  (10)  a  statement  of  the  matrix  modal  equation  of  motion  and  the 
matrix  global  equation  of  motion  at  SEA  are  stated.  Note  the  difference  between  the  loss 
factor  matrix  T\  for  the  modal  and  the  global  equations.  The  difference  lies  merely  in  the 
off-diagonal  elements  and  that  difference  is  accounted  for  by  the  consistency  relationship 
stated  in  Viewgraph  (9).  The  matrices  are  the  transpose  of  each  other. 
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HCS  O  HCS  N-C3 


Modal  SEA 


GJ  E  —  Ke  ;  £  —  { t'a }  ,  —  {Kec 

=  (jla^ccy  ~Vya^~ ^ay))  >  tfa  ~  ^  ^lya 
lobal  SEA 

(0  E  =  ne  ;  E  =  {Ea)  ;  ne  =  \ 

^  rv 


V  =  ’  g  =  {0)5aY) 


Viewgraph  11 


A  recap  is  briefly  stated  in  Viewgraph  (11). 
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The  spatial  dependence  in  the  equation  of  motion 
is  removed  by  introducing  sets  of  orthogonal  modes: 

A  set  of  modes  per  substructure. 

The  phase  dependence  is  removed  by  quadratizing 
linear  quantities: 


Responses  - 

— ► 

Stored  Energy 

Densities 

External  Drives 

— 

External  Input 

Power  Densities 

Dissipations 

- 

Loss  Factors 

Couplings 

Coupling  Loss 

Factors 

Sum  Over  Modes 

Modal  - 

Global 

Modal  Densities 

(Quantity) 

- 

(Am)  (Quantity 

Density) 

V.ll 
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Viewgraph  12 


A  question  may  be  asked  as  to  whether  the  external  input  power  density  into  a 
substructure;  e.g.,  the  (a)  th  substructure,  is  affected  by  the  couplings  to  other 
substructures.  When  SEA  was  developing  under  the  assumption  that  the  coupling 
strengths  are  weak;  e.g., 

flay  Olyy  Vay )  « 1  » 

the  external  input  power  density  (na )  into  the  ( a )  th  substructure  was  assumed  to  be 
unaffected  by  the  couplings;  i.e.,  it  was  assumed  that 

(nea/n°ea)  =  i  • 

Is  this  result  valid  for  coupling  strengths  that  are  strong?  A  heuristic  argument  is 
developed  in  an  attempt  to  answer  this  question. 
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View  graph  13 


For  this  purpose  it  is  found  useful  to  first  define  a  modal  coupling  strength  (^) 

7 

and  its  counterpart,  a  global  coupling  strength  (H«)-  The  definitions  are  cast  in  proper 

y  7 

parametric  forms;  the  modal  and  the  global  coupling  strengths,  (£<£)  and  (Ea), 
respectively,  are  functional  of  only  those  parameters  that  describe  the  substructures  and 
the  couplings  between  them.  These  parametric  forms  are  not  functional  of  the  external 
input  power  into  and/or  the  stored  energies  in  the  substructures.  Also,  note  the  natural 
upper  limits  on  the  modal  and  the  global  coupling  strengths. 
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It  is  convenient  to  define  the  coupling  between  the 
(7)  th  and  the  (a)th  substructures  in  terms  of  a  modal 
coupling  strength.  This  modal  coupling  strength  is 
designated  (££)  and  its  expression  is: 

(ey/e«)  =  Vayiriyy+riay)-1  =  £  <1  . 

Global  coupling  strength 

(Ey/Ea)  =  (nYlna&  =  H«  <  ' 


V.13 
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Viewgraph  14 


The  definitions  of  the  modal  and  the  global  coupling  strengths  are  repeated  and 
the  notion  of  weak  and  strong  couplings  are  explicitly  stated. 
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Modal  coupling  strength 


(er/ea)  =  nay(rir,+narr'  =  £  <1 


Global  coupling  strength 


(Er / Ea)  =  (n7/na)&  =  Era  <  (; nyh 

Weak  coupling: 


e «  1 

E^«  (fly /Ha) 
Strong  coupling: 


i.e.,  r\n  » 77 


lBY<  1 

E£  <(%/««) 


i.e.,  rjn  <  Tj 


lay 


Viewgraph  15 

The  definition  for  the  external  input  power  density,  into  an  isolated  substructure, 
in  terms  of  the  spectral  density  of  the  external  drive  and  the  conductance  is  standard.  The 
expression  for  the  conductance  in  terms  of  the  modal  density  and  the  mass  of  a 
substructure  is  also  standard.  Relying  on  these  standards,  the  conductance  of  a  coupled 
substructure  is  heuristically  suggested.  Of  course,  the  couplings  may  alter  the  modal 
densities  as  well  as  the  masses  that  define  the  substructure.  It  is  assumed  that  such 
alterations  are  second  order  effects  even  when  the  coupling  are  strong;  i.e.,  it  is  assumed 
that  the  na  ’s  and  Ma  ’s  are  robust  parameters  with  respect  to  couplings.  Under  this 
assumption  typically  the  ratio  ( Ga  !G°a)  may  be  determined  as  stated. 
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External  Input  Power 


r^°  c  _  TT° 

^a^fa  LLea  ’ 

G°a  =  (7c/2)(nalMa)  , 

Ga  ^fa  —  . ea ’ 

Ga  =  (n/2)(n'a/Ma) 
nea/n°ea  =  (Ga/G°)  ? 


I  na  =  [na+  X  (nra)]  ;  *4  =  [Ma+L  MyQ 
I  (Ga/G°a)  =  [1+X(S£)][1  +  E  (M./Mjar1 


I 


Viewgraph  16 

i 

Modifications  to  the  global  equation  of  motion  at  SEA  are,  thus,  suggested  so  that 

! 

the  equation  may  be  expressed  in  terms  of  the  external  input  power  density  vector  or 
in  terms  of  the  external  input  power  density  vector  Yl°^-  The  latter  input  power  density  is 
determined  when  the  couplings  are  removed.  How  significant  is  this  modification  to 
SEA?  Before  answering  this  question,  two  simple  examples  in  which  the  modification  is 
insignificant  are  considered  for  the  record.  j 

In  the  first,  weak  coupling  is  assumed.  For  this  case,  weak  coupling  is  defined  as  j 

(X4)  «  land  [I  (MylMa)il\«  1  .  | 

y^a  y^cc 

i 

i 

Under  this  definition  g  =  \ .  \ 
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(Ga/ G°a)  -  [1+S^][1  +Ts(mr/ ma)z/a]  ; 

y*a  y* 2  • 


ra„ 


(GJGZ)  = 


>  1  ;  modal  density -rich  couplings 
1  ;  similar  couplings 
<  1  ;  mass-rich  couplings 


PE  =  (n®)E=ne  =  RU 


(0,/n:)  =  (Ga  /  G°) 


G°asfa  =n°ea  ;  Ga sfa  =nea 


R  =  (RaSar)  ;  /?«=[!  +X ,  Era ]  [1 + yE  (m7  /  ma  0a Y 1 


Y^cc 


[1+  S  (Ml  M)g] 

Y^a  7  a  a 


7i- 1 


(Mr  /  n7 


V.16 
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Viewgraph  17 


In  the  second,  similar  coupling  is  assumed.  In  this  case 
( Ma/na )  =  (My/ tty)  =  (M/n)  =  m 
and,  therefore,  again  R  =  1 .  A  physical  structure  that  trivially  falls  into  this  category 
is  derived  by  introducing  a  “control  line”  on  a  plate;  the  control  line  is  featured  by  the 
dashed  line.  In  this  example  12%  «=  A^_;  «=  A^  and  1.  Therefore,  similar 

couplings  are  assured  and  (Ylecc  /  U°ea  )  =  1.  A  ratio  that  may  be  of  interest  to 
determine  is  (Ea  /  E°a ) .  The  determination  is  presented  in  Viewgraph  (17).  The  results 
indeed  make  physical  sense  and  could  be  apriori  predicted. 
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and  ”>1 

Similar  couplings  and  ( Uea /Yl°ea)-^l 
(Ea  /  Ea )  =  ( Jlaa^t )  ’  V t  =  [^7aa  +  Vyy(Ay 

<V*0  -  (VA>  ^  =  4+Ar 

(^)  l-^iaa-Ax 


n, 


fF  /  F1^  ~  A±diss 

^  in  (a)th  in  [(a)th+(7)th] 


J  ^diss 


If  Tina  ~  riw  >  (Ea  !  K)  =  (At  /  A) 


Viewgraph  18 


What  is  the  relevance  of  the  modified  global  equation  of  motion  at  SEA  to  its 
application,  say,  to  a  noise  control  problem?  The  application  is  illustrated  in 
Viewgraph  18.  In  this  viewgraph,  a  substructure;  e.g.,  the  (a)th  is  depicted  in  isolation 
and  as  coupled  to  another  substructure;  e.g.,  the  (7)  th.  The  coupling  is  introduced  in 
order  to  achieve  a  noise  control  goal,  as  stated  in  Viewgraph  18.  A  generic  definition  of 
a  loss  factor  is  suggested  and  challenged. 
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A  noise  control  goal  is  usually  expressed  in  SEA  by 
a  ratio  of  quantities.  Thus,  a  legitimate  noise  control 
goal  may  be  expressed  in  terms  of  minimizing 

(E«/E‘). 


Isolated 

(Uncoupled) 


Coupled 


Some  of  these  ratios  are  defined  as  loss  factors,  i.e., 

77  =  [ne/(coE)\  ;  What  Ue7 

What  E  ? 


V.18 
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Viewgraph  19 


The  generic  definition  for  a  loss  factor  may  support  a  number  of  definitions.  The 
relationships  among  these  definitions  are  also  presented.  All  these  definitions  pertain  to 
loss  factors  that  are  of  proper  parametric  forms.  Which  of  these  definitions  address  true 
loss  factors  is  also  noted. 

A  loss  factor  requires  accounting  for  the  external  input  power  density  into  a 
structure  (or  a  substructure)  and  the  generated  energy  density  stored  in  the  structure  needs 
to  be  totally  accounted  for.  Otherwise,  questions  such  as  “Where  did  the  energy  go?” 
may  arise  in  a  formalism  that  is  basically  a  statement  of  the  “conservation  of  energy.” 
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Isolated 

(Uncoupled) 


Tlaa=  tn"/(fflE«)  ;  ri“a=rit  =[Uea/((0Ea)]  ; 

%  =  ritY  =  meJ((0Er)]  ;  V"  =  %=  [nea/(©E)] ; 
E  =  (E^+Ey)  ;  Uea  =  ®[(J7aaEa)+(t]r),Er)]  . 

%  =  Ka+VryZin  ;  X  =  TitG&r1  ; 

J7e  =  Vt  [1+SDr1  ;  S£  =  (Ey/E0)  . 


Only  t]aa,  X’  and  %  are  TRUE  loss  factors- 


Viewgraph  20 


A  practical  example  is  presented  in  Viewgraph  20.  A  hollow  rectangular  steel 
structure  is  considered  as  a  master  substructure.  This  substructure  is  coupled  to  another. 
The  second  substructure  is  composed  of  beads  that  fill-in  the  hollow  space.  The  beads 
are  modeled  as  a  fluid  of  low  compressibility  compared  with  that  of  water  and  density 
that  is  comparable  with  that  of  water.  The  low  compressibility  and  the  moderate  density 
yield  a  low  sound  speed;  e.g.,  in  the  range  of  one-fifth  of  the  sound  speed  in  air.  The  low 
speed  of  sound  generates  high  modal  densities.  The  rubbing  of  beads  against  their 
neighbors  yield  high  loss  factors  in  the  beaded-fluid.  The  moderate  density  allows, 
nonetheless,  for  reasonably  high  coupling  strengths.  These  attributes  qualify  the  beaded- 
fluid  substructure  as  a  structural  fuzzy.  As  such,  in  the  steady  state  condition,  a 
significant  portion  of  the  energy  density  is  stored  in  the  adjunct  structure  and  much  of  it 
is  efficiently  dissipated.  In  this  way  the  structural  fuzzy  constitutes  an  effective  damping 
device;  however,  at  best  it  yields  an  overall  loss  factor  that  converges  on  that  of  the 
adjunct  substructure,  as  dictated  by  Le  Chatelier’s  Principle.  Le  Chatelier’s  Principle 
roughly  states:  “When  a  dynamic  system  is  modified  (e.g.,  by  coupling  to  it  a  structural 
fuzzy)  with  the  intention  of  achieving  a  high  degree  of  change  in  the  noise  it  issues,  the 
achievement  is  mitigated  by  the  implementation  of  the  modification,”  e.g.,  the  loss  factor 
converges,  at  best,  to  that  of  the  adjunct  substructure. 
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Structure 


I 


L  =  (10  L,);  L,  =  L2;  h=(L,/20) 


Le  Chatelier’s  principle 


=  <n„/n:„)  mjn,) 

? 

Vt  =  [(iiaa+divyK)]  ;  (nea/iO  =  1 

OWlO  -  (1+3')  [l+Cm/mJElr1 

(■ Ea/E°a )  =  (%M/rie)  [l+(mr / mJElx' 


Viewgraph  21 


The  attributes  of  the  combined  structure,  comprising  the  master  substructure  and 
the  adjunct  substructure,  are  exemplified  in  Viewgraph  21a  and  b.  The  loss  factor  T\t ,  in 
Viewgraph  21c,  is  that  derived  by  setting  the  ratio  {Ylea  /  II °ea )  equal  to  unity. 

[cf.  Viewgraph  20.]  The  convergence  of  the  overall  loss  factor  onto  Tje  is  clearly 
discemable  in  Viewgraph  21. 
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Viewgraph  22 


Variations  on  the  standard  theme  are  depicted  in  this  Viewgraph;  polybeads  and 
beads  in  the  form  of  sand  are  exemplified.  Again,  differentiation  between  7]t  and  r\e  is 
illustrated.  The  convergence  of  the  loss  factor  r\e  on  the  loss  factor  T]yy,  of  the 
structural  fuzzy,  is  of  particular  interest. 
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Viewgraph  23 


In  this  Viewgraph  the  Principle  of  Supplementarity  of  Damping  and  Isolation  is 
annunciated  and  formulated.  This  example  is  the  antonym  of  that  of  Structural  Fuzzy. 
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Another  noise  control  goal  may  be  cast  in  terms  of 
minimizing  the  ratio  of  the  power  density  [(cot]r/)E7] 
dissipated  in  the  passive  (y)th  substructure  to  the 
external  input  power  density  (IIea)  into  the  externally 
driven  (a)th  substructure:  Example  of  the  Principle  of 
Supplementarity  of  Damping  and  Isolation. 


coUri  E  )  +  7?  E  ]  =  IT  ;  E  =  (E  +E  ) 

l\  lyy  yj  ia(x  aJ  ea  ’  v  a  y' 

IJ^CQTlyyEy^  /  Tlecclj  =  ^a^cccc !  Vyy)  a  ~  C EyIEa) 

K  ^  J 

\ 

To  be  minimized 

J 

Increase  Olaa !  Vry)  and  decrease  Z7a 
Principle  of  Supplementarity  of  Damping  and  Isolation 


V.23 
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Viewgraph  24. 


In  this  Viewgraph  the  Principle  of  Supplementary  of  Damping  and  Isolation  is 
exemplified  and  illustrated. 
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